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The  absolute  vorticity  and  angular  momentum  budgets  for 
the  Presidents'  Day  Cyclone,  13-19  February  1979  are  evalu¬ 
ated  employing  the  guasi-Lagrang ian  diagnostic  technique  in 
pressure  coordinates.  The  First  GARP  Global  Experiment 
(FGGE)  Level  Ill-b  data  set  prepared  by  the  European  Center 
for  Medium  Range  Weather  Forecasts  is  used  to  define  the 
storm  and  to  calculate  the  budgets. 

Important  differences  are  found  between  the  FGGE  and 
mass  budget  vertical  motion  estimates.  The  budget  results 
illustrate  the  importance  of  inward  eddy  mode  lateral  trans¬ 
port  and  advection  at  upper  levels  and  mean  mode  transport 
due  to  convergence  at  lower  levels.  The  vertical  redistri¬ 
bution  plays  a  lesser  role  but  is  more  significant  in  the 
angular  momentum  budget.  The  vertical  vorticity  advection 
and  tilting  terms  are  comparable  in  magnitude.  The  largest 
budget  residuals  occur  during  the  period  of  most  intense 
convective  activity  which  indicates  a  need  for  a  better 
understanding  of  the  vertical  redistribution  processes. 
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I. 


INTROD OCTION 


Operational  forecast  models  have  generally  performed 
poorly  in  predicting  the  frequently  observed,  cold-season 
explosive  cyclogenesis  of  maritime  cyclones.  Sanders  and 
Gyakum  (1980)  demonstrated  that  the  National  Meteorological 
Center’s  (NMC)  Limited  Area  Fine-Mesh  Model  (LFM)  predicts 
only  one-fourth  to  one-third  of  the  deepening  rates  of  rap¬ 
idly  deepening  maritime  cyclones. 

The  accurate  forecasting  of  these  meteorological  events 
over  the  open  ocean,  as  well  as  in  coastal  regions,  is  of 
critical  importance  to  both  civilian  and  military  interests. 
The  extreme  conditions  produced  by  these  storms  often  result 
in  less  of  life,  destruction  of  property  and  disruption  of 
operations.  The  recent  sinking  of  the  oil  rig  Ocean  Ranger 
(LeMoyne,  1982)  and  the  attendant  loss  of  life  illustrates 
the  magnitude  of  the  destructive  forces  created  by  intense 
maritime  cyclones.  Throughout  recorded  history,  these 
storms  have  claimed  numerous  ocean-going  vessels  in  prefer¬ 
red  areas  of  explosive  cyclogenesis  such  as  off  the  east 
coast  of  the  0. S.  These  cyclones  often  strike  with  little 
or  no  warning.  More  accurate  and  timely  forecasting  of 
explosively  developing  maritime  cyclones  would  undoubtedly 
reduce  their  impact  on  human  activity. 
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The  failure  of  numerical  forecast  models  in  predicting 
the  intensities  of  vigorous  maritime  cyclones  is  due  to  the 
small  scales  over  which  the  explosive  development  occurs  as 
well  as  the  lack  of  meteorological  observations  over  oceanic 
areas.  This  problem  has  resulted  in  an  inaccurate  represen¬ 
tation  of  the  initial  conditions  in  the  forecast  models. 
The  most  promising  solution  to  this  problem  is  the  incorpo¬ 
ration  of  wind,  temperature,  moisture  and  pressure  data 
obtained  from  satellite  measurements  in  the  operational  pre- 
diction  models.  It  appears  however,  that  a  better  under¬ 
standing  of  the  role  of  dynamic  forcing  in  maritime 
cyclogenesis  is  required  if  this  new  data  source  is  to  be 
fully  utilized. 

This  .study  constitutes  a  portion  of  a  broad-based 
research  effort  to  better  understand  the  processes  leading 
to  the  development,  maturation  and  decay  of  extratropical 
maritime  cyclones.  This  improved  understanding  will  point 
toward  ways  to  improve  the  numerical  prediction  of  maritime 
cyclogenesis. 

More  specifically,  this  thesis  consists  of  circulation 
(vorticity)  and  angular  momentum  budget  analyses  of  the 
Presidents'  Day  Cyclone  (18-19  February  1979);  a  subsynoptic 
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scale  event  which  provided  a  well-defined  example  of 
explosive  maritime  cyclogenesis.  This  storm  has  been  the 
subject  of  several  in-depth  analyses  (Bosart,  1981;  Roman, 
1981;  Uccellini,  Kocin  and  Wash,  1981).  It  produced  record 
breaking  snowfalls  (nearly  60  cm)  in  portions  of  Virginia, 
Maryland  and  Delaware.  The  central  surface  pressure 
decreased  by  more  than  27  mb  during  a  24-hour  period.  The 
NMC  LFM-II  model  under for ecast  the  decrease  in  the  surface 
pressure  minimum  by  an  average  of  8  to  1 6  mb  during  the 
explosive  cyclogenesis  period  (Bosart,  1981).  Thus,  the 
Presidents'  Day  Cyclone  can  be  included  in  a  class  of  explo¬ 
sive  deepeners  not  adequately  handled  by  current  numerical 
forecast  models  (Sanders  and  Gyakum,  1930) .  This  cyclone 
developed  during  one  of  the  First  GARP  Global  Experiment 
(FGGE)  special  observing  periods.  Consequently,  there  are 
additional  observational  data  and  research-based  objective 
analyses  to  aid  in  the  study. 

Roman  (1981)  completed  a  mass  budget  study  of  the  Presi¬ 
dents'  Day  Cyclone  based  on  the  FGGE  data  set.  He  computed 
vertical  velocity  values  in  layers  consistent  with  the  lat¬ 
eral  mass  transport  and  time  tendency  terms.  The  vertical 
velocity  fields  in  the  FGGE  data  set  are  used  to  compute  the 
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vertical  transport  terms  in  the  budget  studies  in  this 
thesis. 

Quasi- La  grangian  diagnostic  (QLD)  procedures  (Wash, 
1978)  are  employed  in  this  study.  The  basic  feature  of  the 
QLD  budget  approach  is  the  computation  of  volume  integrals 
of  atmospheric  properties.  This  volume  is  centered  on  and 
moved  with  the  surface  pressure  minimum  of  the  Presidents' 
Day  Cyclone.  This  approach  incorporates  the  translational 
effects  associated  with  the  cyclone,  as  well  as  the  effects 
due  to  development  or  decay.  The  lateral  exchanges,  sources 
and  sinks,  and  vertical  redistributions  of  meteorological 
variables  are  primarily  the  result  of  developmental  pro¬ 
cesses.  QLD  budgets  yield  quantitative  evaluations  of  the 
exchanges  of  mass,  circulation,  angular  momentum  and  kinetic 
energy  between  the  budget  volume  and  its  environment. 

The  circulation  about  the  surface  cyclone  increases  dra¬ 
matically  during  cyclogenesis,  and  provides  a  convenient 
measure  of  cyclone  intensity.  The  circulation  budget  is 
used  to  evaluate  the  magnitude  of  the  divergence  and  hori¬ 
zontal  and  vertical  advection  terms  in  the  vorticity  equa¬ 
tion.  It  is  well  known  that  upper-level  vorticity  advection 
contributes  to  cyclogenetic  processes  via  the  divergence  and 
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vertical  velocity  fields  (Petterssen,  1955).  This  qualita¬ 
tive  relationship  does  not  however,  provide  answers  to  key 
questions  such  as: 

•  What  is  the  magnitude  and  significance  of  low-level 
horizontal  vorticity  advection? 

•  How  does  the  vertical  advection  of  vorticity  contribute 
to  cyclogenesis? 

•  What  is  the  magnitude  of  the  local  generation  of 
vorticity? 

The  quantitative  analysis  of  circulation  using  budget  diag¬ 
nostics  provides  answers  to  these  questions  through  the 
identification  of  lateral  exchanges,  vertical  redistribu¬ 
tions  and  sources  and  sinks. 

An  alternative,  but  closely  related,  measure  of  vortex 
intensity  is  that  of  angular  momentum.  The  concentration  of 
momentum  about  the  storm  axis  is  primarily  the  result  of 
radial  transport.  Johnson  and  Downey  (1975)  show  that  there 
are  no  sources  of  angular  momentum  within  the  moving  budget 
volume.  Large  scale  forcing  of  storm  angular  momentum 
buildup  occurs  through  inward  transport  of  angular  momentum 
at  lower  levels.  One  objective  of  the  QLD  analysis  of  angu¬ 
lar  momentum  is  to  reveal  the  significance  of  the  vertical 
transport  of  this  property  during  the  development  of  the 
Presidents'  Day  Cyclone. 
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The  budget  analyses  are  done  in  isobaric  vertical  coor¬ 
dinates.  The  choice  of  these  coordinates  is  made  for  con¬ 
venience  in  comparing  the  budget  analyses  with  standard 
diagnostic  and  prognostic  charts.  These  analyses  can  also 
be  converted  to  sigma  vertical  coordinates  over  level  ter¬ 
rain  to  facilitate  comparisons  with  resuLts  from  numerical 
model  predictions. 

In  summary,  the  purpose  of  this  thesis  is  to  examine  a 
well-defined  maritime  cyclogenesis  event  in  terms  of  circu¬ 
lation  and  angular  momentum  from  the  budget  perspective. 
The  specific  objectives  of  this  study  are: 

•  Quantitatively  evaluate  the  horizontal  transport,  ver¬ 
tical  redistribution  and  sources/sinks  of  vorticity  and 
angular  momentum  within  a  budget  volume  centered  on  the 
Presidents'  Day  Cyclone; 

•  Compare  the  vertical  velocity  structure  of  the  Presi¬ 
dents'  Day  Cyclone  derived  from  the  mass  budget  verti¬ 
cal  transport  (Roman,  1981)  with  that  of  the  F3GE  data 
set . 

A  brief  description  of  the  FGGE  data  set  and  a  synoptic 
overview  of  the  Presidents'  Day  Cyclone  are  presented  in  the 
following  chapter.  Chapter  III  contains  an  abbreviated 
description  of  the  QLD  technique.  A  more  detailed  explana¬ 
tion  of  the  quasi- La grangian  coordinate  system  and  general¬ 
ized  budget  equations  is  presented  in  Appendix  B.  The 
results  of  the  budget  analyses  are  presented  in  Chapter  IV. 
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Chapter  V 
mendations 


contains  the  conclusions  of  this  study  and  recom- 
for  further  research. 
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II.  DATA  BASE  AND  SYNOPTIC  OVERVIEW 

A.  DATA  BASE 

The  synoptic  descriptions  and  budget  analyses  of  the 
Presidents*  Day  Cyclone  are  based  solely  on  the  FGGE  Level 
Ill-b  data  set.  More  specifically,  FGGE  analyses  between 
1200  GMT  18  February  and  1200  GMT  20  February  1979  are  used 
to  define  the  development  and  movement  of  this  storm.  The 
FGGE  analyses  were  prepared  by  the  European  Center  for 
Medium  Range  Weather  Forecasts  ( ECMWF)  using  a  three-dimen¬ 
sional  optimal  interpolation  scheme  on  a  high  resolution 
(1.875  degrees  of  latitude  and  longitude)  nonstaggered  grid. 
This  nonlinear  analysis  scheme  produced  high-guality,  verti¬ 
cally  and  horizontally  consistent  mass  and  momentum  fields. 
A  more  detailed  description  of  the  analysis  scheme  is  found 
in  the  Global  Weather  Experiment .  The  Daily  Global  Analysis, 
Part  J,  (1981),  published  by  ECMWF. 

The  validity  of  the  FGGE  data  set  in  defining  the  Presi¬ 
dents'  Day  Cyclone  was  demonstrated  by  Roman  (1981) .  His 
subjective  comparisons  with  the  Bosart  (1981)  analyses,  the 
Fleet  Numerical  Oceanography  Center  analyses  and  the  NMC  LFM 
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analyses  clearly  illustrate  the  accuracy  of  the  FGGE  data  in 
representing  the  atmospheric  conditions  daring  the  period. 

Some  of  the  advantages  in  using  the  FGGE  data  set  as 
opposed  to  conventional  analyses  are: 

♦  The  incorporation  of  a  variety  of  data  sources  includ¬ 
ing  off-tame  and  late  data; 

•  The  high  resolution  regular  grid  is  particularly  useful 
in  examining  the  subsynoptac  seals  features  of  the 
Presidents'  Day  Cyclone; 


•  The  data  set  includes  several  derived  fields  such  as 
relative  humidity  and  vertical  velocity; 

•  The  analyses  fully  encompass  the  area  of  development 
and  movement  of  the  Presidents'  Day  Cyclone. 

The  FGGE  Level  Ill-b  data  set  includes  fields  of  pres¬ 
sure  height,  u  and  v  wind  components,  temperature  and  verti¬ 
cal  velocity  for  all  standard  pressure  levels  and  relative 
humidity  for  all  standard  levels  up  to  300mb  (Table  I). 
Additionally,  absolute  vorticity  and  thickness  fields  were 
derived  for  use  in  this  thesis. 

B.  SYNOPTIC  OVERVIEW 

The  purpose  of  this  section  is  to  review  the  changes  in 
the  mass  and  momentum  fields  over  the  eastern  U.S.  and  west¬ 
ern  Atlantic  that  led  to  the  development  of  the  Presidents' 
Day  Cyclone.  This  review  differs  from  previous  discussions 
of  the  storm  (Bosart,  1981;  and  Uccellini,  Kocin  and  Wash, 
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1981)  by  focusing  on  these  changes  in  relation  to  the  budget 
volume  translating  with  the  surface  pressure  minimum.  The 
discussion  is  based  on  the  pressure,  temperature  and  wind 
analyses  and  the  derived  vorticity  and  thickness  fields  from 
the  FGGE  data  set.  The  discussion  in  terms  of  vorticity  and 
thermal  advection  serves  to  provide  a  background  for  the 
circulation  and  angular  momentum  budget  study  analyses  pre¬ 
sented  in  Chapter  IV. 

All  of  the  synoptic  figures  in  this  thesis  depict  budget 
volumes  with  lateral  boundaries  at  six  and  nine  degrees 
latitude.  The  six-degree  boundary  was  selected  because  it 
generally  includes  the  intense  part  of  the  storm.  The 
nine-degree  boundary  was  chosen  since  it  encompasses  nearly 
all  of  the  peripheral  features  related  to  the  storm.  The 
six  and  nine  degree  circles  are  referred  to  throughout  this 
section  as  the  inner  and  outer  boundaries.  The  circular 
representation  of  the  budget  volume  lateral  boundaries 
introduces  a  small  distortion  since  these  boundaries  should 
be  slightly  elliptical  on  the  conformal  projections  used  in 
the  figures.  Also,  the  position  of  the  budget  volume  at 
1200  GMT  18  February  is  based  on  the  first  distinguishable 
surface  pressure  minimum  of  the  Presidents'  Day  Cyclone  (not 
shown) . 
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TABLE  I 

FGGE  LEVEL  III-b  DATA  FIELDS  FOR  TIMES  SHOWN  BELOW 


SURFACE  PRESSURE  FIELDS  ARE  ALSO  AVAILABLE 


LEVEL  (mb  x  10) 


100  85  70  50  40 


30  25  20 


15  10 


PRESSURE  HEIGHTS 
VERTICAL  VELOCITY 
RELATIVE  HUMIDITY 
TEMPERATURE 
U  WIND  COMPONENTS 
V  WIND  COMPONENTS 


X  X  X  X  X 
X  X  X  X  X 
X  X  X  X  X 
X  X  X  X  X 
X  X  X  X  X 
X  X  X  X  X 


X  X  X  X  X 
X  X  X  X  X 
X 

X  X  X  X  X 
X  X  X  X  X 
X  X  X  X  X 


TIME  PERIODS  AVAILABLE  (GMT) 


m  r» — 17  n  g — n — So — & — t2 

18  FEB  19  FEB  20  FEB 


Figure  1  depicts  the  6-hourly  positions  of  the  storm's 
surface  pressure  minima  prior  to,  during  and  following  the 
period  of  explosive  cyclogenesis  from  0000  GMT  19  February 
through  0000  GMT  20  February  1979.  During  this  time,  the 
surface  pressure  minima  decreased  by  27  mb  for  an  average 
deepening  rate  in  excess  of  1  mb  per  hour  for  24  hours.  The 
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maximum  6-hour  decrease  was  19.3  mb  between  1200  GMT  and 


1800  GMT  19  February.  Figs.  2  through  6  depict  the  surface 
pressure  patterns  and  the  850  and  500mb  height  contours  and 
isotherms  at  12-hour  intervals  from  1200  GMT  18  February 
through  1200  GMT  20  February.  The  lOOOmb  thermal  fields 
appear  with  the  surface  pressure  patterns  since  the  surface 
temperature  fields  are  not  included  in  the  FGGE  Level-IIIb 
data  set.  The  500mb  wind  vectors  and  absolute  vorticity 
fields  at  12-hour  intervals  are  presented  in  Figs.  7  and  8. 
The  examination  of  these  figures  yields  qualitative  esti¬ 
mates  of  mid-level  vorticity  advection.  Finally,  the 
1000-500mb  thickness  patterns  and  lOOOmb  height  contours 
appear  in  Figs.  9  and  10.  These  figures  are  useful  in 
determining  lower  tropospheric  thermal  advection. 

At  1200  GMT  18  February  the  dominant  surface  synoptic 
feature  is  a  massive  cold  high  pressure  system  centered  over 
northern  New  York  with  a  ridge  to  the  south  over  and  east  of 
the  Appalachians  (Fig.  2a).  This  ridge  is  a  hydrostatic 
reflection  of  the  shallow  cold  air  entrenched  east  of  the 
Appalachians.  The  thermal  trough  at  the  lOOOmb  level  in  this 
region  is  located  within  the  northwestern  portion  of  the 
outer  budget  volume  boundary.  Another  feature  of  interest 
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Figure  1.  Presidents  Day  Cyclone  track  at  6-hour  intervals 
from  1800  GMT  18  February  through  1200  GMT  20 
February,  1979  with  5  degree  radius  circles. 

Times  are  date  and  hour  in  GMT  (Z) .  Central 
surface  pressures  in  mb  appear  below  times. 


23 


at  the  surface  is  the  inverted  pressure  trough  in  the 
central  Gulf  of  Mexico.  Additionally,  a  *eak  trough  is  seen 
near  the  center  of  the  budget  volume.  The  850  mb  contour 
pattern  (Fig.  2b)  shows  a  strong  ridge  over  the  northeastern 
states  that  extends  into  the  Great  Lakes  region,  and  a  sharp 
trough  from  North  Dakota  southeastward  across  the  inner 
budget  volume  boundary.  The  850mb  isotherm  pattern  is  pri¬ 
marily  zonal  in  the  vicinity  of  the  budget  volume  with  lit¬ 
tle  indication  of  the  underlying  shallow  cold  air.  The  key 
feature  at  the  500mb  level  (Fig.  2c)  is  a  short  wave  embed¬ 
ded  in  a  nearly  zonal  flow  across  the  upper  Midwest.  The 
wave  is  generally  in  phase  with  the  thermal  trough  at  this 
level.  A  weak  ridge  is  seen  in  the  50Dmb  contour  pattern 
north  of  the  inner  volume  boundary.  The  absolute  vorticity 
maximum  (Fig.  7a)  extends  into  the  northern  portion  of  the 
budget  volume  and  southward  into  Texas.  Positive  vorticity 
advection  (PVA)  is  occurring  to  the  west  and  northwest  of 
the  budget  volume  from  Illinois  southwest  ward  into  Louisi¬ 
ana.  The  1000-500mb  thickness  pattern  (Fig.  9a)  indicates  a 
broad  area  of  strong  cold  advection  throughout  the  lower 
troposphere  over  the  western  North  Atlantic.  Weak  cold 
advection  is  present  across  the  northeastern  budget  volume 
bounda  ry . 
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Figure  2.  a.  Surface  pressure  and  lOQOmb  isotherms,  b. 

850mb  height  contours  and  isotherms  and  c.  500mb 
contours  and  isotherms  for  1200  GMT  18  February 
1979.  Surface  pressure  is  p-IOOOmb.  Contours  are 
labeled  in  decameters  and  isotherms  (°C)  are 
dashed.  Circles  are  6  and  9  degree  budget  volume 
lateral  boundaries. 
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Significant  changes  occur  in  the  low-level  patterns 
during  the  12  hour  period  ending  at  0000  GMT  19  February. 
The  surface  high  pressure  center  (Fig.  3a)  moves  southeast¬ 
ward  to  Massachusetts.  More  importantly,  however,  is  the 
development  of  a  cyclonic  circulation  off  the  coast  of  South 
Carolina.  There  was  little  indication  12  hours  earlier  that 
cyclogenesis  would  commence  in  this  area.  One  might  argue 
that  the  inverted  trough  previously  i-n  the  Gulf  could  have 
signaled  the  onset  of  cyclogenesis,  if  the  trough  had  moved 
to  the  current  position  of  the  low  with  the  500mb  flow. 
However,  the  1800  GMT  surface  analysis  (not  shown)  shows  no 
indication  of  the  trough  translating  northeastward  across 
the  Florida  peninsula.  Coinciding  with  the  cyclonic  devel¬ 
opment  off  the  coast  is  a  reorientation  and  strengthening  of 
the  lOOOmb  thermal  gradient  along  the  coast  from  Georgia  to 
Maryland  within  the  inner  volume  boundary.  The  presence  of 
the  coastal  front  and  its  contribution  to  the  cyclone  devel¬ 
opment  is  discussed  by  Bosart  (1981).  In  addition,  a  sig¬ 
nificant  thermal  ridge  has  developed  within  the  inner  volume 
boundary. 

The  primary  change  in  the  850mb  contour  pattern  (Fig. 
3b)  is  the  southeastward  extension  of  the  trough  over  the 
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surface  low.  This  trough  now  encompasses  the  entire  north¬ 
west  quadrant  of  the  budget  volume.  Little  change  is  seen 
in  the  850mb  isotherm  pattern  during  the  12-hour  period. 
The  thermal  patterns  at  the  1300  and  350mb  levels  indicate 
low-level  destabilization  within  the  inner  volume  boundary 
which  would  favor  organized  convection. 

The  500mb  wave  (Fig.  3c)  has  moved  eastward  to  Illinois 
and  is  more  in  phase  with  the  thermal  trough  at  this  time. 
The  wave  is  beginning  to  penetrate  the  northwest  quadrant  of 
the  cuter  budget  volume  boundary.  A  weak  ridge  is  seen  in 
the  isotherm  pattern  at  this  level  throughout  the  inner  vol¬ 
ume  boundary.  The  vorticity  pattern  (Fig.  7b)  shows  weak 
P7A  over  the  northern  half  of  the  volume.  Thermal  advec¬ 
tion,  as* depicted  in  the  thickness  pattern  (Fig.  9b),  is 
confined  to  the  inner  volume  radius  with  warm  advection  to 
the  east  and  weak  cold  advection  to  the  west  of  the  cyclone 
center. 

The  1200  GMT  19  February  surface  pressure  pattern  (Fig. 
4a)  indicates  significant  cyclogenesis  has  occurred.  The 
minimum  pressure  decreased  by  9  mb  in  che  12-hour  period 
ending  at  this  time.  The  surface  isobars  show  a  highly  sym¬ 
metrical  circular  cyclone  structure  located  entirely  within 
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the  inner  budget  radius.  The  1003mb  thermal  ridge  is  con¬ 
siderably  sharper  and  is  in  nearly  perfect  alignment  with 
the  cyclone  center  and  associated  pressure  trough.  The 
lOOOmb  isotherm  pattern  reveals  the  development  of  a  baro- 
clinic  zone  to  the  southwest  of  the  cyclone  within  the  inner 
budget  volume.  The  850mb  height  contours  (Fig.  4b)  now  show 
closed  cyclonic  circulation  centered  in  the  northwest  quad¬ 
rant  of  the  inner  volume  radius  which  indicates  a  northwest 
tilt  of  the  cyclone  axis  toward  the  cold  air.  Significant 
changes  in  the  850mb  isotherm  pattern  have  occurred  during 
this  12-hour  period.  The  development  of  a  thermal  ridge 
over  the  surface  center  indicates  that  the  warm  core  is  pen¬ 
etrating  to  higher  levels.  This  warming  is  likely  the 
result  of  latent  heat  release  through  convection  and  the 
warm  advection  shown  at  the  850mb  level  12  hours  earlier. 

The  500mb  wave  (Fig.  4c)  is  now  centered  over  Pennsylva¬ 
nia  and  is  nearly  in  phase  with  the  thermal  trough  within 
the  northwest  quadrant  of  the  outer  budget  volume  boundary. 
The  500mb  vorticity  center  (Fig.  7c)  is  now  within  the  outer 
radius  with  significant  PVA  in  the  inner  radius  primarily  to 
the  northeast  of  the  surface  center.  The  thickness  pattern 
(Fig.  9c)  shows  a  warm  ridge  centered  on  the  cyclone  with 
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warm  adv ection  throughout  the  eastern  half  and  cold  advec- 
tion  throughout  the  western  half  of  the  outer  budget  radius. 

Intense  cyclogenesis  continues  during  the  next  12-hour 
period  ending  at  0000  GMT  20  ?ebruary.  During  this  period 
the  surface  pressure  minimum  (Fig.  5a)  decreases  by  nearly 
18  mb.  At  the  same  time,  the  cyclone  takes  a  more  easterly 
course  following  the  maximum  500mb  PVA  region.  The  vertical 
structure  indicates  that  the  cyclone  is  nearing  the  mature 
stage.  The  circular  surface  and  850mb  centers  are  nearly 
coincident  and  the  isotherm  patterns  show  a  vertical  align¬ 
ment  of  the  thermal  ridge  through  the  850mb  level  (Fig.  5b). 
The  500mb  wave  (Fig.  5c)  is  almost  directly  over  the  surface 
cyclone  and  is  encompassed  by  the  outer  budget  radius.  The 
thermal  trough  now  lags  the  500mb  wave  with  an  intensified 
thermal  ridge  within  the  eastern  half  of  the  outer  volume 
boundary.  Thus,  the  warm  core  now  extends  through  the  low 
troposphere.  The  thickness  pattern  (Fig.  10a)  clearly 
illustrates  a  mature  warm  core  structure.  The  vertical 
extension  of  the  warm  core  through  the  530mb  level  is  most 
likely  due  to  convection  since  little  thermal  advection  is 
seen  at  this  level  prior  to  0003  GMT.  The  mature  structure 
indicates  that  little  further  deepening  of  the  cyclone  can 
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be  expected.  The  500mb  vorticity  pattern  (Fig.  8a)  shows 
PV A  throughout  the  eastern  half  of  the  budget  volume  with 
maximum  PVA  due  east  of  the  cyclone  center. 

The  surface  pressure  pattern  (Fig.  6a)  at  1200  GMT  20 
February  shows  that  the  cyclone  center  has  accelerated  in  an 
easterly  direction  during  the  previous  12  hours  with  little 
further  deepening.  The  storm  center  moves  nearly  250  nm 
eastward  for  an  average  velocity  of  approximately  20  knots. 
Little  change  is  seen  in  the  thermal  structure  through  the 
low  troposphere  in  the  inner  budget  volume.  This  might 
indicate  that  a  moist  adiabatic  lapse  rate  has  been  achieved 
by  convection  throughout  much  of  the  inner  budget  radius. 
The  500mb  vorticity  pattern  (Fig.  9b)  now  shows  moderate  to 
strong  P7A  in  the  eastern  half  and  definite  negative  vortic¬ 
ity  advection  (NVA)  in  the  northwest  quadrant  of  the  outer 
volume  boundary. 

The  intent  of  this  synoptic  overview  has  been  to  present 
the  highlights  of  the  development  of  the  Presidents*  Day 
Cyclone  in  terms  of  the  mass  and  circulation  fields  as  a 
prelude  to  the  storm-following  budget  analyses  in  this  the¬ 
sis.  The  goal  of  budget  studies  of  this  storm  is  to  quan¬ 
tify  the  forcing  mechanisms  of  explosive  maritime 
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cyclogenesis.  A  brief  discussion  of  the  quasi-Lagra ng ian 
technique  employed  in  the  budget  study  is  presented  in  the 
next  chapter. 
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Figure  7 


500m b absolute  vorticity  fields  (x  10-s  Sec-M 
fS4oW19d  Jfctors  for  a*  1200  GMT  18  February 
1^9,  bj  0000  GMT  19  February  1979  and  c.  1$00 
GMT  i 9  February  1979.  Circles  are  6  and  9  degree 
budget  volume  lateral  boundariss. 
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Figure  9. 


1 000-500mb  thickness  and  1 OOQmb  contour  Datt erns 
for  a.  1200  GMT  18  February  1979.  b.  0000  GMT  19 
February  1979  and  c.  1200  GMT  19  February  1979. 
Thickness  contours  are  dashed  and  labeled  in 
decameters.  1000  mb  contours  are  labeled  in 
meters.  Circles  are  5  and  9  degree  budcret  volume 
lateral  boundaries. 
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III.  QUASI-LAGRANGIAN  DIAGNOSTIC  TECHNIQUE 

The  objective  of  budget  studies  is  to  evaluate  the 
changes  of  atmospheric  properties  within  a  specific  atmos¬ 
pheric  volume.  Normally  these  studies  have  been  used  prima¬ 
rily  to  study  the  general  circulation  and  planetary  scale 
balances  of  atmospheric  properties  such  as  mass,  momentum 
and  kinetic  energy.  For  synoptic-scale  studies  the  approach 
has  been  to  evaluate  changes  in  these  properties  within  a 
stationary  volume  (Rung  and  Baker,  1975;  Smith,  1973).  The 
effects  of  synoptic-scale  disturbances  migrating  through  the 
volume  are  quantitatively  evaluated  relative  to  the  fixed 
volume  only.  This  approach  is  not  well  suited  to  the  study 
of  cyclone  development  since  changes  in  the  budget  volume 
can  be  due  to  the  movement  of  the  cyclonic  circulation  into 
or  out  of  the  volume  as  well  as  due  to  cyclone  development. 

The  QLD  budget  approach  employed  in  this  thesis  uses  a 
budget  volume  which  translates  with  the  cyclone.  This 
method  permits  a  quantitative  evaluation  of  the  changes  in 
the  budget  property  relative  to  the  cyclone.  Specifically, 
this  system  isolates  the  horizontal  advection  and  transport 
relative  to  the  cyclone  from  the  horizontal  advection  and 
transport  associated  with  the  translation  of  the  cyclone. 
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Further,  the  vertical  redistribution  of  the  budget  property 
is  the  result  of  the  vertical  transport  processes  in  the 
vicinity  of  the  cyclone.  Thus,  the  QLD  technique  is  well 
suited  to  the  study  of  cyclogenesis.  Table  II  presents  an 
updated  list  of  recently  completed  studies  employing  the  QLD 
method. 


TABLE  II 

PREVIOUS  QU ASI-LAG RANG IAN  DIAGNOSTIC  BUDGET  STUDIES 
(EXTRACTED  IN  PART  FROM  MASH,  1978) 


Property 

Available  Potential  Energy 

Circulation 

Mass 

Absolute  Angular  Momentum 

Kinetic  Energy 

Absolute  Angular  Momentum 

Mass 

Mass 


Researcher  (s) 

Spaete  1974 
Wash  1975 

Johnson  and  Downey  1976 
Johnson  and  Downey  1976 
Chen  and  Bosart  1977 
Hash  1978 
Roman  1981 
Tallman  1982 


The  budget  equation  is  evaluated  with  respect  to  a 
budget  volume  in  spherical  coordinates  (Appendix  B) .  In 
this  thesis,  the  budget  volume  is  centered  over  and  main¬ 
tains  its  position  relative  to  the  surface  pressure  minimum 
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of  the  Presidents*  Day  Cyclone.  The  volume  is  defined  in 
the  vertical  by  10  mandatory  pressure  levels  (1000,  850, 

700,  500,  400,  300,  250,  200,  150  and  100  mbs)  and  horizon¬ 
tally  by  radii  for  example,  six  degrees  of  latitude.  Appen¬ 
dix  B  establishes  the  quasi-Lagrangian  budget  volume 
coordinate  system  and  the  generalized  budget  equations  that 
can  be  applied  to  any  selected  atmospheric  property. 

The  first  step  in  the  calculation  of  the  budgets  is  to 
interpolate  the  u  and  v  wind  component  values  at  the  rectan¬ 
gular  gridpoints  in  the  FGGE  data  set  to  the  budget  volume 
gridpoints.  The  interpolated  winds  are  then  converted  to 
normal  and  tangential  components  relative  to  the  circular 
budget  volume  boundaries.  Additionally,  the  movement  of  the 
budget  volume  (H0)  is  used  to  compute  the  normal  and  tangen¬ 
tial  components  relative  to  the  translating  volume. 

The  budget  equations  for  the  circulation  and  angular 
momentum  budgets  in  pressure  coordinates  are  presented  in 
Tables  III  and  IV.  In  the  vorticity  budget  formulation,  the 
horizontal  advection  and  divergence  terms  are  combined  to 
form  the  lateral  transport  using  the  Gauss  Divergence 


The  vertical  redistribution  occurs  through  the  vertical 
advection  term  and  the  friction  and  tilting  terms  complete 
the  budget.  The  angular  momentum  budget  is  formulated  by 
combining  the  continuity  and  angular  momentum  equations. 
Since  there  is  no  divergence  term,  the  vertical  redistribu¬ 
tion  occurs  through  a  vertical  transport  term.  The  friction 
term  completes  this  budget. 

The  budget  volume  grid  consists  of  36  equally  spaced 
points  around  each  circular  lateral  boundary  at  each  pres¬ 
sure  level.  Absolute  vorticity  and  angular  momentum  values 
are  computed  at  each  of  the  volume  gridpoints.  The  calcula¬ 
tions  of  the  time  tendency,  transport  and  source/sink  terms 
are  based  on  these  gridpoint  values.  The  tilting  term  in 
the  vorticity  equation  is  computed  as  a  source.  There  are 
no  sources  of  absolute  angular  momentum  within  the  budget 
volume,  and  the  only  sink  is  due  to  friction  which  is  com¬ 
puted  using  a  stability  independent  parameterization  scheme 
(Johnson  and  Downey,  1976)  with  a  constant  geostrophic  drag 
coefficient  (0.035). 

The  transport  terms  are  broken  into  mean  and  eddy  modes 
(Table  V).  The  purpose  of  this  statistical  partitioning 
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is  to  isolate  changes  in  the  budget  property  associ¬ 
ated  with  the  synoptic-scale  mean  and  eddy  circulations. 
The  mean  and  eddy  modes  represent  the  transport  by  the  irro- 
tational  and  rotational  wind  components,  respectively.  In 
this  context,  the  mean  and  eddy  modes  are  analogous  to  the 
divergence  and  horizontal  advection  terms  in  the  vorticity 
eguation.  In  this  study  the  divergence  and  advection  terms 
are  calculated  for  comparison  with  the  mean  and  eddy  modes. 
The  comparisons  serve  to  further  substantiate  the  validity 
of  the  mean  and  eddy  mode  partitioning  in  the  QLD  technique. 

The  FGGE  vertical  velocity  fields  were  calculated  based 
on  divergences  derived  through  a  nonlinear  normal  mode  ini¬ 
tialization  scheme.  An  objective  of  this  thesis  is  to  com¬ 
pare  the  FGGE  vertical  velocities  with  the  vertical  velocity 
values  derived  via  the  QLD  method  in  the  mass  budget  of  the 
Presidents*  Day  Cyclone  (Roman,  1981).  The  purpose  of  this 
comparison  is  to  aid  in  the  assessment  "of  the  vertical 
velocities  computed  in  the  budget  technique.  The  comparison 
is  presented  in  conjunction  with  the  analysis  of  the  budget 
vertical  transport  term  in  the  next  chapter. 
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TABLE  III 


QU  ASI-  LAGR  ANGIAN  CIRCULATION  BUDGET  EQUATION  IN  ISOBARIC 

VERTICAL  COORDINATES 

(after  Mash,  1978) 


lOOmb  B  2it  i 

C  »  f  f  f  — C  r2sinB  dadB(-dp) 

a  lOOOmb  00  9  a 


where  C  is  absolute  circulation  and  £  is 
a  a 

absolute  vorticity. 


The  budget  equation  is 


~a=LT(C  )  +VT ( C  )+S(C  ) 

o  t  cl  cl  » 


where  the  lateral  transport  is 
lOOmb  2ir  i 

LTCC  )=  /  /  -(U-W)CrsinB  do(-dp). 

a  lOOOmb  0  g  B  a  B 


and  the  vertical  advection  is 


VTCC  1=  r2sinB  dadB. 

a  0  0  9  3P 


The  source/sink  term  is 


lOOmb  B  2ir  i  dc  ■> 

SCC  )  =*  I  /  /  —  ~a  r2sin8  dodB(-dp). 

a  lOOOmb  00  g  at 


See  Appendix  B  for  notation  explanation. 
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TABLE  IV 


QUAS I- LAGRANGI AN  ANGULAR  MOMENTUM  BUDGET  EQUATION  IN 
y  ISOBARIC  VERTICAL  COORDINATES 


(after  Wash,  1978) 


lOOmb  S  2j  i  „ 

G  «  I  f  f  — a  r2 sinB  dadS(-dp) 

a  lOOOmb  00  g  a 


where  G  is  the  volume  integral  of  absolute  angular 

cl  ^ 

momentum  g  . 

a 


The  budget  equation  is 
5  Q 

a*»LT  ( G  )  +VT  (G)+S(G) 
ot  a  a  a 


where  the  lateral  transport  is 

lOOmb  2n  1  • 

LT (G  )-  /  /  -<&-») Rgrsin$  do(-dp) | 

a  lOOOmb  0  g 


and  the  vertical  divergence  of  the  vertical  transport  is 
lOOmb  3  2  it  15 

VT(G  )«  /  /  /  — r-(ug  )r2sin8  dadB  C-dp). 

a  lOOOmb  00  g3P  a 


The  sink  term  is 


S(G  ) 
a 


lOOmb  8 
/  / 
lOOOmb  0 


2  IT 

/ 

0 


I  *£a 

g  dt 


r2sin8 


dadB  (-dp). 


See  Appendix  3  for  notation  explanation. 
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TABLE  V 

PARTITION  OF  LATERAL  TRANSPORT  INTO  MEAN  AND  EDDY  MODES 


(after  Wash  ,  1978) 


lOOmb  2ir 


C*(U-w)*  rsinB  da(-dp)|fl 

Si  p  p 


EM  ( C  )  =  /  / 

^  1  Art A mK  A 


lOOOmb  0 


lOOmb  2ir  i  . .  .-a  .  , 

MM ( C  ) =  /  /  —  C  rsinB  doC-dp) |ft 

a  lOOOmb  0  ^  a  0  p 


lOOmb  2ir  i  „  . 

EM ( G  )=  /  /  -  g'(U-W)*  rsinB  da(-dp) 

a  lOOOmb  0  9  a  6  8 


lOOmb  2 it  i  _a  - -  . 

MM  ( G  )  **  /  /  —  g=  CU-W)  fl  rsinB  daC-dp)|. 

a  lOOOmb  0  9  a  6  “ 


EM  is  the  eddy  mode  and'  MM  is  the  mean  mode. 


Note  that 


and  (  )  (  )  -  (  ) a  is  the  deviation  of 


the  property  from  its  mean  around  a  lateral 
boundary.  See  Tables  III  and  IV  and  Appendix 
B  for  notation  explanation. 
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IV.  SOASI^LAGR ANGIAN  BUDGET  ANALYSIS 

A.  GENERAL 

The  absolute  circulation  or  vorticity  and  absolute  angu¬ 
lar  momentum  budget  results  are  presented  in  this  chapter. 
Emphasis  is  placed  on  the  vorticity  budget  since  the  results 
are  similar  and  vorticity  analyses  are  more  commonly  used  to 
describe  cyclone  development.  The  angular  momentum  results 
are  discussed  in  depth  where  significant  differences  exist, 
particularly  in  the  vertical  transport  term. 

The  results  are  presented  by  means  of  layer  or  level/ 
time  or  time  period  cross-sections  and  selected  vertical 
profiles.  Results  are  presented  for  two  radii,  6°  and  9°, 
which  were  illustrated  in  the  synoptic  discussion.  Values 
are  displayed  on  levels  or  layers  over  6-hour  time  periods. 
For  example,  the  600mb  layer  extends  from  700  to  500mb  and 
the  1500  GMT  time  period  runs  from  1200  to  1800  GMT.  The 
sign  convention  is  positive  for  processes  which  produce  a 
vorticity  or  angular  momentum  increase  such  as  inward  trans¬ 
port,  inward  horizontal  advection  and  sources.  Upward  ver¬ 
tical  velocity  (omega)  is  negative  (pressure  coordinates). 
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B.  TIME  TENDENCY 


The  dramatic  buildup  of  absolute  vorticity  within  the 
budget  volume  during  the  development  of  the  Presidents'  Day 
Cyclone  is  illustrated  by  the  vorticity  time  sections  in 
Fig.  11.  These  patterns  reflect  the  vorticity  maximum 
developing  within  the  budget  radii  at  the  1000  mb  level 
(Fig.  12)  and  the  vorticity  maxima  moving  into  the  budget 
circles  at  the  500  and  300  mb  levels  (Figs.  7  and  13 
respectively) . 

The  vorticity  time  section  patterns  (Fig.  11)  within 
both  radii  are  similar  with  minima  occurring  in  the  225mb 
layer  during  the  1500  GMT  18  February  time  period.  These 
minima  are  seen  in  terms  of  anticyclonic  shear  and  curvature 
within  the  budget  radii  in  the  300mb  wind  field  (Fig.  14a). 
The  vorticity  values  are  greater  in  the  6°  radius  than  those 
in  the  9°  radius  after  the  0900  GMT  19  February  period.  The 
larger  gradient  with  respect  to  time  in  the  inner  radius 
depicts  the  rapid  increase  and  concentration  of  cyclonic 
circulation  as  the  storm  develops.  Maxima  in  both  volumes 
occur  in  the  upper  layers  during  the  0900  GMT  20  February 
period.  The  larger  vorticity  increase  with  time  in  the 
upper  layers  is  due  to  the  development  of  cyclonic  shear  and 
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curvature  as  the  upper  tropospheric  trough  and  associated 
jet  maximum  dominate  the  southern  half  of  the  inner  budget 
volume  (Fig.  14b).  The  vorticity  profiles  (Fig.  15)  at  the 
beginning  and  end  of  the  intense  cyclogenesis  period  show 
that  the  largest  increase  occurs  in  the  upper  layers. 

Vorticity  time  tendencies  computed  over  6-hour  intervals 
are  depicted  in  Fig.  16.  These  patterns  show  positive  ten¬ 
dencies  in  all  layers  and  time  periods  with  the  exception  of 
the  lower  layers  within  the  inner  radius  during  the  final 
time  period.  Both  patterns  show  primary  and  secondary  max¬ 
ima.  The  primary  maximum  in  the  9°  radius  (350mb  layer, 
0300  GMT  19  February  period)  occurs  earlier  than  the  corre¬ 
sponding  maximum  in  the  6°  radius.  The  outer  region  of  the 
cyclone  is  spinning-up  earlier  in  response  to  the  influence 
of  the  upper  level  trough  and  jet  maximum  at  this  time.  As 
the  trough  overtakes  the  lower  level  cyclone,  the  6®  radius 
shows  a  large  vorticity  increase.  The  secondary  maxima 
occur  during  the  same  period  (2100  GMT  19  February)  in  each 
volume,  but  at  a  higher  layer.  These  maxima  are  associated 
with  the  strong  cyclonic  shear  and  curvature  north  of  the 
upper  level  jet  core  which  is  located  within  the  budget  vol¬ 
ume  at  this  time  (Fig.  14b).  The  magnitude  of  the  secondary 
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9°  VOLUME 


6°  VOLUME 


Fig  are  11. 


18FEB  19FEB  20  FEB 


Absolute  vorticity  time  sections  (x  10-5 
sec-*)  . 
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Figure  12.  lOOOmb  wind  vectors  and  absolute  vorticity  (x 

10”s  sec-1)  for  a.  3000  GMT  19  February  1979.  b. 
1200  GMT  19  February  1979  and  c.  0000  GMT  20 
February  1979.  Circles  are  6  and  9  degree 
budget  volume  lateral  boundaries. 
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Figure  13. 


300  mb  wind  vectors  and  and  absolute  vorticit 
(x  10~5  sec-*r  for  a.  0600  GMT  19  February  an 
b.  1800  GMT  19  February  1979.  Circles  are  6 
9  degree  budget  volume  lateral  boundaries. 


I 

and 
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Figure  14. 


300  mb  wind 
a.  1800  GMT 
February  1979. 

budget  volume  lateral  boundaries 


vectors  and 
18  February 
Circles 


lsotachs 
and  b.  06 
are  6  and 


sec-1 1 
GMT  20 
degr 
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Figure  15. 


absolute  vorticity  profiles  for  radius  six 
(solid)  and  radius  nine  (dashed)  . 


maximum  in  the  6°  radius  is  nearly  twice  that  in  the  9° 
radius,  because  the  maximum  cyclonic  curvature  and  shear 
vorticity  occur  within  the  inner  boundary. 

The  absolute  angular  momentum  time  section  patterns 
(Fig.  17)  are  similar  to  those  of  absolute  vorticity,  except 
the  values  in  the  outer  radius  exceed  those  in  the  inner 
radius.  This  difference  is  because  storm  angular  momentum, 
%  =JS  •  l  *  (0-j^,  is  proportional  to  the  distance  (R)  from 
the  storm  center.  The  absolute  angular  momentum  (Fig.  18) 
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9°  VOLUME 


6°  VOLUME 


Figure  16. 


18FEB  1 9  FEB  20  FEB 


Absolute  vorticity  time  tendencies  over  6  hours 
(x  1 0~ 10  sec-2).  Dashed  indicates  negative 
values. 
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and  vorticity  time  tendency  patterns  exhibit  remarkable  sim- 


ilarities  with  virtually  coincident 

maxima. 

The 

larger 

angular 

momentum  tendencies  in  the 

larger 

volume 

again 

ref lect 

the  dependence  of  this  prope 

rty  on 

radius . 

The 

buildup  of  angular  momentum  within  the  budget  volume  por¬ 
trayed  in  these  patterns  is  consistent  with  the  findings  of 
Johnson  and  Downey  (1976). 

In  summary,  the  vorticity  and  angular  momentum  patterns 
show  that  the  largest  increases  occur  within  the  24-hour 
period  of  intense  cyclogenesis  (0000  GMT  19  February  through 
0000  GMT  20  February).  The  maximum  vorticity  buildup  occurs 
in  the  upper  layers  because  the  polar  jet  core  moves  into 
the  budget  volume.  The  primary  difference  between  the  vor¬ 
ticity  and  angular  momentum  buildup  is  that  the  maximum 
angular  momentum  is  found  at  a  lower  layer  particularly  in 
the  6°  radius  volume  (Fig.  17). 

C.  LATERAL  TRANSPORT 

The  lateral  transport  of  absolute  vorticity  (Fig.  19)  is 
the  most  important  forcing  term  in  the  budget.  The  maximum 
inward  transport  occurs  in  the  upper  layers,  although  the 
transport  is  positive  for  nearly  all  of  the  layers.  The 
maximum  in  the  smaller  volume  occurs  later,  is  in  a  lower 
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layer  and  is  considerably  greater  in  magnitude  than  the  max¬ 
imum  in  the  larger  volume.  Both  maxima  occur  within  the 

24-hour  period  of  intense  cyclogenesis;  in  particular,  the 
larger  maximum  in  the  inner  volume  occurs  during  the  period 
of  most  rapid  deepening  (1200-1800  GMT  19  February).  The 
maximum  transport  in  the  larger  volume  occurs  earlier  (0600 
GMT  19  February)  as  the  vorticity  maximum  approaches  the 
northwestern  portion  of  the  outer  boundary  (Fig.  13a) . 
Twelve  hours  later,  the  vorticity  maximum  crosses  the  6° 
boundary  (Fig.  13b)  and  produces  the  larger  maximum  trans¬ 
port  at  the  inner  radius. 

As  described  in  the  preceding  chapter,  the  total  lateral 
transport  can  be  partitioned  into  mean  and  eddy  mode  compo¬ 
nents.  The  purpose  of  this  partitioning  is  to  break  t-he 
total  transport  into  components  to  improve  the  interpreta¬ 
tion  of  the  lateral  transport  dynamic  forcing.  The  mean 
mode  represents  the  mean  vorticity  transport  by  the  average 
normal  wind  component  (mean  divergence  or  convergence) , 
while  the  eddy  mode  is  the  transport  due  to  the  covariance 
of  the  vorticity  and  normal  wind  component  deviations.  In 
more  general  terms,  the  mean  mode  transport  is  due  to  the 
mean  inflow  and  outflow  of  the  cyclone  and  the  eddy  mode 
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9°  VOLUME 


6°  VOLUME 


Figure  17. 


18FEB  19FEB  23  FEB 


Absolute  angular  momentum  time  sections  (x  10l° 
cm  sec  -2)  . 
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9°  VOLUME 


6°  VOLUME 


Figure  18. 


18FEB  19FEB  20  FEB 


Absolute  angular  momentum  time  tendencies  (x  10s 
cm2  sec-2) .  Dashed  indicates  negative  values. 
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9°  VOLUME 


6°  VOLUME 


Figure  19. 


18FEB  19FEB  2DFEB 


Absolute  vorticity  lateral  transport  (x  10-l° 
see-2) .  Dashed  lines  indicate  negative  values. 
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transport  reflects  transport  due  to  the  asymmetries  in  the 
cyclone  and  jet  streak. 

An  alternative  partitioning  of  the  lateral  transport 
term  follows  the  traditional  vorticity  equation  with  the 
horizontal  advection  and  divergence  terms.  The  area-aver¬ 
aged  absolute  vorticity  horizontal  advection  was  calculated 
for  comparison  with  the  eddy  mode  transport.  The  divergence 
term  was  not  calculated  in  this  study  and  its  comparison 
with  the  mean  mode  transport  must  be  accomplished 
indirectly. 

The  vorticity  eddy  mode  transport  time  sections  (Fig. 
20)  are  very  similar  to  the  total  transport,  except  the  mag¬ 
nitudes  are  larger.  The  maximum  in  the  smaller  volume  again 
occurs  during  a  later  time  period  and  is  greater  in  magni¬ 
tude  than  the  outer  radius  maximum.  The  primary  features  in 
the  area-averaged  horizontal  vorticity  advection  patterns 
(Fig.  21)  generally  correspond  to  those  in  the  eddy  mode 
patterns.  The  upper  level  maxima  in  the  inner  volume  occur 
during  the  same  time  period,  although  the  P V A  maximum  occurs 
in  a  higher  layer.  The  large  maximum  in  the  upper  tropo¬ 
sphere  results  from  the  movement  of  the  strong  vorticity 
maximum  and  jet  streak  into  the  cyclone  region  as  described 
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in  the  synoptic  discussion  (Figs.  7,  8,  13  and  14).  The 
similarities  in  these  patterns  suggest  a  strong  correlation 
between  the  eddy  mode  and  the  traditional  view  of  the  role 
of  PVA  at  upper  levels  in  the  development  of  extratropical 
cyclones  (Petterssen,  1956). 

The  most  significant  difference  between  the  eddy  mode 
and  area-averaged  advection  patterns,  particularly  within 
the  inner  volume,  is  the  PVA  (approximately  15  x  10~10 
sec-2)  in  the  925  mb  layer  during  the  2100  GMT  19  February 
period.  This  low-level  PVA  does  not  agree  with  the  advec¬ 
tion  implied  by  the  lOOOmb  wind  vectors  and  absolute  vortic- 
ity  patterns  (Fig.  12).  These  patterns  actually  indicate  a 
net  NVA  within  the  inner  radius  since  the  wind  vectors  are 
parallel  or  cross  the  vorticity  isopleths  toward  higher  val¬ 
ues.  In  an  effort  to  resolve  this  discrepancy,  the  absolute 
vorticity  advection  values  at  each  gridpoint  were  examined. 
The  data  indicate  that  the  low-level  PVA  is  occurring  prima¬ 
rily  at  the  inner  radii  (less  than  6°)  to  the  east  and 
northeast  of  the  storm  center.  The  low-level  PVA  in  the 
time  sections  is  the  result  of  asymmetries  in  the  low-level 
vorticity  contours  and  winds  relative  to  the  translating 
storm  volume.  The  effect  of  the  translating  volume  is  to 
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enhance  the  low-level  advecting  (relative)  winds  (easter¬ 
lies)  and  to  reduce  the  upper  level  advecting  winds  (wester¬ 
lies)  .  This  accounts  for  the  smaller  difference  in  magni¬ 
tude  between  the  upper  and  lower  level  vorticity  advection 
maxima  in  a  guasi-Lagrang ian  compared  to  a  Eulerian  budget 
study.  Aside  from  these  small  discrepancies,  the  overall 
similarities  in  the  patterns  support  the  close  relation 
between  the  eddy  mode  lateral  transport  and  horizontal 
advection. 

The  mean  mode  vorticity  patterns  (Fig.  22)  show  outward 
lateral  transport  in  the  mid  to  upper  troposphere  and  inward 
transport  in  the  lower  layers.  These  patterns  reflect  vor¬ 
ticity  transport  due  to  the  low-level  convergent  and  upper- 
level  divergent  cyclone  structure.  The  upper-level  outflow 
in  the  smaller  volume  increases  with  time  to  the  end  of  the 
period.  In  contrast,  the  smaller  maximum  outflow  value  for 
the  larger  volume  is  achieved  relatively  early  and  is  main¬ 
tained  thereafter.  The  larger  mean  mode  inflow  within  the 
6°  radius  volume  is  the  result  of  larger  vorticity  values 
within  the  inner  boundary  as  well  as  the  more  vigorous  mass 
circulation  at  the  inner  radii. 
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9°  VOLUME 


6°  VOLUME 


Figure  20. 


18FEB  19FEB  2D  FEB 


Absolute  vorticity  eddy  mode  lateral  transport 
(x  10” 10  sec-2) . 
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9°  VOLUME 


6°  VOLUME 


Figure  2 1 . 


18FEB  19FEB  2D  FEB 


Areal  average  absolute  vorticity  horizontal 
advection  (x  10~l°  sec-2)  . 
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9°  VOLUME 


60  VOLUME 


Figure  22. 


18FEB  19FEB  23  FEB 


Absolute  vorticity  mean  mode  lateral  transpor 
(x  10~l°  sec-2).  Dashed  lines  indicate  negat 
values. 


67 


The  vorticity  lateral  transport  profiles  (Fig.  23) 
illustrate  the  mean  and  eddy  mode  contributions  to  the  total 
lateral  transport.  The  mean  mode  outflow  (corresponding  to 
divergence)  is  overwhelmed  by  the  eddy  mole  inflow  (PVA)  in 
the  upper  layers  while,  in  the  lower  layers  the  mean  mode 
inflow  (convergence)  dominates.  The  relationship  between 
the  mean  mode  transport  and  the  divergence  is  also  illus¬ 
trated  by  noting  that  the  level  of  the  zero  value  in  the 
mean  mode  profile  (690mb)  agrees  with  the  mid-tropospheric 
level  of  nondivergence  (LND)  implied  by  the  lateral  mass 
transport  patterns  (Roman,  1981)  . 

The  absolute  angular  momentum  lateral  transport  patterns 
(Figs.  24-26)  are  also  similar  to  the  corresponding  vortic¬ 
ity  patterns.  Here  again,  the  larger  transport  values  in 
the  outer  radius  are  the  result  of  the  storm  angular  momen¬ 
tum  values  increasing  with  the  radius.  The  primary  differ¬ 
ence  within  the  inner  volume  is  that  there  is  outward 
transport  in  the  mid  to  upper  layers  during  the  final  time 
period  in  Fig.  24b  whereas  the  transport  is  inward  in  Fig. 
19b.  This  indicates  that  the  mean  mode  contribution  to  the 
lateral  transport  is  greater  in  the  angular  momentum  budget 
than  in  the  vorticity  budget. 
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Figure  23.  Total  (solid)  ,  mean  (dashed)  and  eddy  (dash-dot) 
absolute  vorticity  lateral  transport  profiles 
for  the  0000-0600  3HT  19  February  period. 

In  summary,  the  lateral  transport  term  is  the  most 
important  term  in  the  vorticity  and  angular  momentum  budg¬ 
ets.  The  comparisons  between  the  eddy  and  mean  mode  parti¬ 
tions  and  the  horizontal  advection  show  that  the  cyclone's 
mean  lower  level  inflow  (convergence)  forces  the  low  tropo¬ 
sphere  vorticity  and  angular  momentum  buildup.  At  the  upper 
levels,  the  eddy  mode  inward  transport  or  horizontal  advec¬ 
tion  is  predominant  and  is  responsible  for  the  upper  tropos¬ 
pheric  increase  of  vorticity  and  angular  momentum  within  the 
volume. 
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9°  VOLUME 


6°  VOLUME 


Figure  24. 


18FEB  19FEB  20FEB 


Absolute  angular  moaentum  lateral  transport  (x 
10s  cm2  sec-2)  .  Dashed  lines  indicate  negative 
values. 
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9°  VOLUME 


6°  VOLUME 


Figure  25. 


18FEB  1 9  FEB  20FEB 


Absolute  angular  momentum  mean  mode  lateral 
transport  (x  10s  cm2  sec-2).  Dashed  lines 
indicate  negative  values. 
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9°  VOLUME 


6°  VOLUME 


Figure  26. 


18FEB  1 9 FEB  20  FEB 


Absolute  angular  momentum  eddy  mode  lateral 
transport  (x  10s  cm2  sec”2) . 
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D.  VERTICAL  REDISTRIBUTION 

The  vertical  redistribution  terms  (vertical  vorticity 
advection  and  the  divergence  of  the  vertical  angular  momen¬ 
tum  transport)  were  calculated  using  the  vertical  velocity 
fields  in  the  FGGE  data  set  (Fig.  27).  The  method  used  to 
calculate  these  values  is  briefly  described  in  the  preceding 
chapter.  Roman  (1981)  calculated  vertical  velocities  con¬ 
sistent  with  the  lateral  transport  and  time  tendency  terms 
in  the  QLD  mass  budget  of  the  Presidents*  Day  Cyclone  (Fig. 
28)  ;  hereafter  referred  to  as  the  QLD  vertical  velocities. 
Note  that  the  QLD  vertical  velocities  are  mean  values  for 
each  layer  in  the  budget  volume  and  can  only  estimate  a  mean 
mode  vertical  redistribution  if  used  in  the  budget  calcula¬ 
tions.  The  time  sections  of  the  FGGE  and  QLD  vertical 
velocities  reveal  significant  differences.  The  FGGE  values 
are  generally  greater  with  maximum  vertical  motions  that 
occur  earlier  and  at  lower  levels  (600mb  vice  400mb) . 

The  FGGE  vertical  velocity  values  during  the  period  of 
most  rapid  deepening  (Fig.  29)  exceed  the  QLD  values,  espe¬ 
cially  below  400mb.  The  discrepancy  between  the  FGGE  and 
mass  budget  vertical  motion  profiles  is  surprising.  If  one 
assumes  that  the  FGGE  vertical  motions  ace  derived  from  the 
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9°  VOLUME 


6°  VOLUME 


Figure  28. 
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OLD  vertical  velocity  time  sections  (Mbac 
sec-1) .  Dashed  linas  indicate  negative  values. 
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adiabatic  version  of  the  nonlinear  normal  mode  initializa¬ 


tion,  the  mass  budget  omegas  might  be  expected  to  be  larger 
since  they  include  diabatic  contributions.  More  study  is 
required  to  discover  the  causes  for  these  differences  in 
vertical  motion  estimates. 

The  vertical  vorticity  advection  time  sections  (Fig.  30) 
reveal  intricate  but  similar  patterns  at  each  radius.  Max¬ 
ima  occur  nearly  at  the  level  of  the  maxima  in  the  FGGE  ver¬ 
tical  velocity  fields  (used  to  calculate  this  term)  but 
during  a  later  time  period.  Weak  downward  advection  is 
indicated  in  the  mid-layers  with  weak  upward  advection  in 
the  upper  layers.  The  mid-layer  downward  advection  must  be 
attributed  to  a  reversal  in  the  vertical  vorticity  gradient 
since  no  downward  tropospheric  vertical  velocities  are  seen 
in  the  layer/area  averaged  FG3E  data  and  a  weak  vorticity 
decrease  with  height  is  evident  in  the  low  to  mid-layers  in 
the  radius  six  vorticity  time  sections  (Fig.  11).  The  maxi¬ 
mum  value  of  the  vorticity  vertical  advection  term  (7  x 
10-*°  sec-2)  is  smaller  by  a  factor  of  two  than  the  time 
tendency  term  (13  x  10“10  see”2)  .  These  terms  however,  are 
considerably  smaller  than  the  maximum  value  in  the  total 
lateral  transport  term  (25  x  10~10  sec-2)  .  Thus,  the 
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contribution  by  the  vertical  advection  term  to  the  budget 
results  is  relatively  small. 

In  the  angular  momentum  budget,  the  vertical  redistribu¬ 
tion  occurs  through  the  vertical  divergence  of  the  vertical 
angular  momentum  transport  (Fig.  31) .  Significant  differ¬ 
ences  exist  between  the  inner  and  outer  volume  patterns. 
Within  radius  six,  a  net  angular  momentum  loss  is  indicated 
in  the  two  lowest  layers,  particularly  after  the  0300  GMT  19 
February  time  period.  However,  within  the  9°  volume,  the 
pattern  reveals  a  gain  in  the  lowest  layers  and  a  loss  in 
the  450mb  layer  during  the  most  intense  period  of  surface 
development  (1200-1800  GMT  19  February).  Both  volumes  have 
mid- layer  gains  throughout  the  48-hour  period,  with  signifi¬ 
cant  gains  in  the  upper  layers  in  the  inner  volume  that 
reach  a  maximum  during  the  final  period.  The  larger  volume 
pattern  shows  a  net  loss  in  the  upper  layers  throughout  the 
period  of  study  with  a  maximum  during  the  1500  GMT  19  Febru¬ 
ary  period. 

The  angular  momentum  vertical  transport  (Fig.  32) 
through  each  level  was  calculated  to  aid  in  resolving  the 
primary  features  of  the  vertical  divergance  term.  These 
patterns  show  that  the  maximum  upward  transport  during  the 
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Figure  29. 


FGGE  (solid) 
profiles  for 


and  QLD  (dashed) 
1200-1300  GMT  19 


vertical  velocity 
February  1979. 
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9°  VOLUME 


6°  VOLUME 


Figure  30. 
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Vorticity  vertical  advection  (x  10-10  sec-2). 
Dashed  indicates  negative  or  downward  advection. 
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initial  time  period  rises  from  850mb  to  the  700mb  level  by 
the  0900  GMT  19  February  period  in  the  6°  radius.  More 
importantly  however,  is  the  dramatic  height  increase  of  the 
zero  isopleth  during  the  early  stages  of  cyclogenesis  in  the 
inner  radius.  The  depth  of  the  upward  angular  momentum  ver¬ 
tical  transport  reaches  a  peak  during  the  most  intense 
development  period.  Thus,  a  vertical  divergence  (loss)  of 
angular  momentum  occurs  in  the  low  troposphere  while  a  con¬ 
vergence  (gain)  of  this  property  is  found  in  the  mid  and 
upper  troposphere. 

The  9°  pattern  is  considerably  more  detailed.  The  most 
striking  feature  is  the  downward  transport  through  all  lev¬ 
els  during  the  1500  and  2100  GMT  19  February  periods.  This 
indicates  the  presence  of  downward  vertical  velocities  (not 
seen  in  the  area-average  omegas)  in  the  outer  radius  donut. 
These  downward  values  are  related  to  the  split  in  the  maxi¬ 
mum  of  area-averaged  FGGE  vertical  velocities  (Fig.  27)  and 
the  smaller  vertical  transport  values  in  the  outer  radius 
compared  to  the  inner  radius.  The  product  of  the  downward 
velocities  and  the  larger  angular  momentum  values  at  outer 
radii  is  large  enough  to  overcome  the  6°  radius  area-aver¬ 
aged  upward  transport  and  produce  the  downward  transport  in 
the  9°  pattern  during  this  interval. 
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In  a  general  sense,  the  relative  contribution  of  the 
vertical  redistribution  term  (vertical  divergence  of  the 
vertical  transport)  to  the  angular  momentum  budget  can  be 
seen  by  comparing  maximum  values.  Phe  largest  value  of  this 
term  in  the  inner  volume  (12  x  10s  cm2  see-2)  is  two  thirds 
of  the  maximum  inner  volume  time  tendency  term  (18  x  10s  cm2 
sec-2).  As  in  the  vorticity  budget,  the  lateral  transport 
is  the  dominant  term  with  a  maximum  value  of  35  x  10s  cm2 
sec-2  in  the  inner  volume. 

The  analyses  in  this  section  have  shown  the  upward  ver¬ 
tical  redistribution  of  vorticity  and  angular  momentum  from 
the  low  to  the  upper  troposphere.  In  the  vorticity  budget, 
the  large  lower-tropospheric  vorticities  are  advected  upward 
although  this  term  plays  a  relatively  small  role  in  the 
budget.  In  the  angular  momentum  budget,  the  vertical  diver¬ 
gence  of  the  transport  is  generally  negative  in  the  lower 
layers  and  positive  aloft  and  plays  a  more  important  role  in 
the  budget.  The  FGGE  vertical  motion  estimates  were  found 
to  differ  significantly  from  the  mass  budget  estimates. 
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Divergence  of  absolute  angular  momentum  vertical 
transport  (x  10s  cm2  sec-2) .  Dashed  lines 
indicate  negative  values  (losses). 
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Figure  32.  Absolute  angular  momentum  vertical  transport  (x 
107  mb  cm2  sec-2) .  Dashed  indicates  negative 
values  (upward  transport). 
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E.  SOURCES/SINKS 


The  vorticity  and  angular  momentum  sinks  due  to  fric¬ 
tional  dissipation  were  computed  in  the  lowest  layer  using  a 
stability  independent  parameterization  scheme  (Johnson  and 
Downey,  1976)  to  ascertain  the  significance  of  this  contri¬ 
bution  to  the  budgets.  The  calculated  values  were  found  to 
be  negligible  compared  to  the  other  budget  terms  and  thus 
have  no  significant  effect  on  the  budget  residuals. 

The  vorticity  tilting  term  (Fig.  33)  was  calculated  and 
its  values  are  generally  comparable  in  magnitude  to  the  time 
tendency  and  vertical  advection  values.  The  patterns  show 
weak  cyclonic  vorticity  generation  in  all  layers  during  the 
early  time  periods  and  continuing  in  the  mid-layers  through 
the  final  period.  Anticyclonic  vorticity  generation  is 
indicated  in  the  upper  and  lower  layers  during  the  1500  and 
2100  GMT  19  February  time  periods.  Strong  pattern  similari¬ 
ties  exist  between  the  tilting  and  vertical  advection  (Fig. 
30)  terms  during  the  later  stages  of  cyclogenesis.  During 
this  time,  these  terms  exhibit  maxima  of  opposite  sign  in 
the  lower  layers  and  sign  reversals  in  the  mid-layers.  An 
analysis  of  the  components  of  the  tilting  term  might  reveal 
the  reasons  for  the  similarities  and  differences  between  the 
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tilting  and  vertical  advection  patterns.  In  this  study  how¬ 
ever,  the  tilting  term  was  calculated  only  to  determine  its 
significance  relative  to  the  primary  budget  terms.  In  any 
case,  the  tilting  term  patterns  are  consistent  with  the 
cyclone  study  of  Chen  and  Bosart  (1979) .  Overall,  the  pat¬ 
terns  reveal  a  net  positive  contribution  during  the  early 
stages  of  cyclogenesis  and  a  net  negative  contribution  dur¬ 
ing  the  latter  stages. 

F.  RESIDUAL 

Residuals  were  computed  to  determine  the  integrity  of 
the  budgets.  The  vorticity  budget  residuals  (Fig.  34)  are 
near  zero  in  the  lower  layers  within  the  5°  volume  while  in 
the  9°  volume  the  values  are  relatively  small  throughout. 
The  large  negative  values  (a  vorticity  deficiency)  in  the 
upper  layers  in  the  inner  volume  during  the  6-hour  period  of 
explosive  cyclogenesis  (1200  to  1800  GMT  19  Feb)  deserve 
more  intense  scrutiny.  The  extreme  value  in  the  225mb  layer 
(approximately  -15  x  10~l°  sec-2)  is  nearly  60*  of  the  maxi¬ 
mum  in  the  total  lateral  transport  term  (nearly  25  x  10-l° 
sec  -2)  in  the  inner  radius.  An  equally  large  positive  res¬ 
idual  (vorticity  excess)  is  found  in  the  highest  layer  dur¬ 
ing  the  0300  GMT  20  February  period.  Errors  introduced  by 
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interpolation,  finite  differencing,  and  spatial  and  time 
averaging  can  contribute  to  the  residuals.  These  contribu¬ 
tions  however,  should  be  reasonably  consistent  during  all 
time  periods  and  in  all  layers. 

Since  the  largest  negative  residuals  occur  during  the 
6-hour  period  of  most  explosive  cyclogenesis  and  deep  con¬ 
vection  (Bosart,  1981),  it  seems  likely  that  the  intense 
mesoscale  convective  processes  are  not  adequately  resolved 
in  the  FGGE  data  set.  As  in  most  synoptic  studies,  a  major 
uncertainty  lies  in  the  vertical  velocity  fields.  The  pres¬ 
ence  of  large  residuals  suggests  that  the  normal  mode  verti¬ 
cal  velocities  (FGGE)  do  not  explain  the  vertical  vorticity 
redistribution  during  the  period  of  explosive  deepening. 

The  angular  momentum  budget  residual  patterns  (Fig.  35) 
are  more  complex  than  the  vorticity  residual  patterns.  The 
major  imbalances  are  again  found  during  the  period  of 
intense  development  and  convective  activity,  1200  to  1800 
GMT  19  February,  in  the  outer  volume.  The  large  negative 
values  indicate  an  angular  momentum  deficit  in  the  mid  and 
upper  layers.  Within  the  6°  volume,  a  positive  maximum 
(excess)  is  found  in  the  350mb  layer  during  the  0300  GMT  20 
February  period.  However,  the  largest  inner  radius  deficits 
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Tilting  term  in  the  vorticity  budget  (x  10~l° 
sec-2).  Dashed  lines  indicate  negative  values. 
Contour  interval  is  0.8  units. 
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are  found  in  the  lowest  layers  and  reach  a  maximum  during 
the  final  time  period.  This  indicates  a  problem  in  the 
parameterization  of  the  frictional  torque  for  this  data  set. 

An  interesting  point  common  to  both  budgets  is  that  the 
maximum  upper  level  deficits  occur  in  the  time  interval 
between  the  double  maxima  in  the  tendency  patterns  (Figs.  16 
and  18).  At  first  glance  this  might  indicate  a  phase  prob¬ 
lem  in  the  budget  calculations.  Closer  examination  however, 
reveals  that  a  doubling  of  the  time  tendency  in  the  interval 
between  the  twin  maxima  is  required  to  compensate  for  the 
vorticity  and  angular  momentum  deficits  in  the  upper  levels. 
Furthermore,  this  compensation  would  require  a  substantial 
increase  in  the  magnitudes  of  the  existing  maxima  in  the 
tendency  patterns.  Thus,  it  is  unlikely  that  a  problem  in 
the  budget  calculations  is  the  primary  cause  of  this  incon¬ 
sistency.  Again,  the  vertical  motion  estimates  (FGGE)  used 
in  the  budget  computations  appear  suspect,  particularly 
since  the  largest  deficits  common  to  both  budgets  occur  dur¬ 
ing  the  period  of  most  rapid  deepening. 
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Figure  34.  Vorticity  budget  residuals  (x  10“l°  sec 
Dashed  lines  indicate  negative  values, 
interval  is  3.0  units. 


89 


N  CJ 


18FEB 


19FEB 


20  FEB 


9°  VOLUME 


6°  VOLUME 


Figure  35.  Angular  momentum  budget  residuals  (x. 10s 
sec-2)  .  Dashed  lines  indicate  negative 
Contour  interval  is  8.0  units. 
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V. 


CONCLUSIONS  AND  RECOMMENDATIONS 

The  absolute  vorticity  and  angular  momentum  budgets  for 
the  Presidents*  Day  Cyclone,  18-19  February  1979  were  exam¬ 
ined  employing  the  quasi-Lagrang ian  diagnostic  technique  in 
pressure  coordinates.  The  FS3E  Level  Ill-b  data  set  pre¬ 
pared  by  the  ECMWF  was  used  to  define  the  development  and 
movement  of  the  storm  as  well  as  to  calculate  the  budgets. 

The  budget  calculations  in  this  study  have  yielded  a 
number  of  interesting  results.  The  lateral  transport  is  the 
most  significant  term  in  each  budget.  Maximum  values  of  the 
lateral  transport  are  generally  twice  those  of  the  time  ten¬ 
dencies,  which  is  usually  the  second  largest  term.  The  mag¬ 
nitudes  of  the  vertical  redistribution  terms  are 
approximately  half  those  of  the  time  tendency.  The  fric¬ 
tional  dissipation  was  small  while  the  vorticity  tilting 
term  is  equal  in  magnitude  to  the  vorticity  vertical  advec- 
tion.  The  residuals  in  each  budget  were  largest  during  the 
period  of  most  intense  convective  activity  and  explosive 
cyclogenesis  within  the  larger  volume  in  the  angular  momen¬ 
tum  budget  and  within  the  smaller  volume  in  the  vorticity 
budget. 
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The  mass  budget  average  vertical  motion  estimates  were 
compared  with  the  area-averaged  vertical  motion  from  the 
ECMWF  analyses.  Significant  differences  were  found.  Sur¬ 
prisingly,  the  ECMWF  vertical  motion  estimates,  which  are 
believed  to  be  from  an  adiabatic  version  of  the  nonlinear 


normal  mode  initialization,  were  larger  than  the  essentially 
kinematic  average  vertical  motions  from  the  mass  budget. 

The  mean  and  eddy  mode  partitioning  of  the  lateral 
transport  and  the  correlations  with  the  horizontal  advection 
and  divergence  components  of  the  lateral  transport  show  that 
the  mean  low-level  inflow  supplies  vorticity  and  angular 


momentum  to  the  low-troposphere.  However,  the  upper  level 
horizontal  advection  or  eddy  mode  inward  transport  is  prima¬ 
rily  responsible  for  the  absolute  vorticity  and  storm  angu¬ 
lar  momentum  buildup  in  the  upper  troposphere  during  the 
development  of  the  Presidents'  Day  Cyclone. 

Several  recommendations  for  further  research  emerged 
during  the  development  of  this  thesis: 


•  Vertical  velocity  fields  should  be  computed  based  on 
the  FGGE  Level  Ill-b  data  set  using  the  kinematic 
method  in  an  attempt  to  incorporate  tne  diabatic  heat¬ 
ing  contribution.  Calculation  of  the  budgets  based  on 
these  vertical  velocity  fields  could  then  be  accom¬ 
plished  in  an  effort  to  reduce  the  residuals  and  more 
accurately  evaluate  the  budget  terms; 


•  Vertical  cross-sections  of  the  budget  terms  through  the 
storm  center  would  allow  a  more  detailed  analysis  of 
the  large-scale  dynamic  forcing  mechanisms; 
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•  The  divergence  term  of  the  vorticity  budget  .should  be 
analyzed  to  further  understand  the  relationship  bexween 
the  advective  and  transport  forms  of  the  lateral  trans¬ 
port; 


•  Energy  diagnostics  (kinetic  and  available  potential) 
should  also  be  completed  on  this  cyclone  for  a  more 
comprehensive  description  of  the  explosive  cyclogenesis 
exhibited  by  this  class  of  cyclones. 
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APPENDIX  A 


•  cm 

•  ECMWF 

•  FGGE 

•  GARP 

•  GMT 

•  LFM 

•  m 

•  mb 

•  NVA 

•  NMC 

•  PVA 

•  QLD 

•  sec 


S  AND  ACRONYMS 


COMMON  ABBREVIATION 

centimeters 

European  Center  for 
Forecasts 

First  GARP  Global  E 
Global  Atmospheric 
Greenwich  Mean  Time 
Limited  area  Fine-M 
meters 
millibars 

Negative  Vorticity 
National  Meteorolog 
Positive  Vorticity 
Quasi -Lagrangian  Di 
seconds 


Medium  Range  Weather 

xper iment 
Research  Program 


esh 

Advection 
ical  Center 
Advection 
agnostics 
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APPENDIX  B 


QUASI-LAGRANGIAN  BUDGET  FRAMEWORK  AND 
GENERALIZED  BUDGET  EQUATIONS 


The  purpose  of  this  appendix  is  to  present  the  formula¬ 
tion  of  the  storm  budget  volume  coordinate  system  and  the 
generalized  budget  equations  which  can  be  applied  to  any 
desired  atmospheric  property.  The  information  contained 
herein  is  extracted  from  Wash  (1 978)  as  adopted  from  Johnson 
and  Downey  (1975A). 

Figure  36  depicts  the  storm  budget  volume  in  spherical 
coordinates.  The  following  notation  is  relevant  to  this 
figure: 

a  Azimuthal  coordinate 


{$  Angular  radial  coordinate 

r  Position  vector  from  earth's 

~  center 


R  Position  vector  from  storm  center 

~  to  any  point  in  budget  volume 


m,l,k  Orthogonal  unit  vectors  along 

~  ~  ~  meridional,  azimuthal  and 

vertical  directions  of  spherical 
coordinates  for  storm 


U 


Wind  velocity  relative  to  earth 
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w 

~  o 

Horizontal  velocity  of  reference 
axis 

w 

Horizontal  velocity  at  any  point 
on  the  volume  boundary. 

Table  VI 

contains  a  listing  of  the  generalized  budget 

equations.  The  following  is  an  explanation  of  the  notation 
in  these  equations: 


(U-ff) 

~ 

Horizontal  wind  component  relative 
to  budget  volume  that  is  normal  to 
lateral  boundary 

<jJ 

Vertical  velocity  in  pressure 
coordinates 

r  sin  p 

Radius  of  budget  volume 

P 

Pressure 

g 

gravit  y 
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ri-C 


An. 


Figure  36.  Quasi-Langrangian  storm  budget  volume  coordinate 
system  (extracted  from  Mash,  1978). 
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TABLE  VI 


GENERALIZED  QUASI- LA GR AN GIAN  BUDGET  EQUATION  IN  PRESSURE 

COORDINATES 


(after  Hash,  1978) 


lOOmb  B  2ir  - 

p»  /  /  /  —  fr2sin&  dadS(-dp) 

lOOOmb  0  0  g 

where  F  is  the  volume  integral  of  the  desired  budget 
property  f. 


The  budget  equation  is 

=  LTCF)  +VT(F)  +S  (F) 


St 


where  the  lateral  transport  is 


LT(F) 


LT(F)  ■( JU-W)  g  f rsinB  da(-dp)| 

B 

and  the  vertical  redistribution  is 


lOOmb  B 

VT(Fl»  /  /  / 

lOOOmb  0  0 


g  3p 


The  source/sink  term  is 


lOOmb  B 

S(F)=*  /  /  f 

lOOOmb  0  0 
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